Electron heat diffusivity in radially-bounded ergodic region of toroidal plasma

Kanno, Ryutaro*; Nunami, Masanori*; Satake, Shinsuke*; Matsuoka, Seikichi; Takamaru, Hisanori*

Nuclear Fusion, 58(1), p.016033_1 - 016033_7, 2018/01

https://doi.org/10.1088/1741-4326/aa94ed

The electron heat transport in a torus plasma which involves a radially-bounded ergodic region, where flux surfaces are partially destroyed by perturbative magnetic fields, is studied. In this paper, we have demonstrated that the radial heat conduction by the particles' parallel motion is reduced by trapped particles.

Benchmark of the local drift-kinetic models for neoclassical transport simulation in helical plasmas

Huang, B.*; Satake, Shinsuke*; Kanno, Ryutaro*; Sugama, Hideo*; Matsuoka, Seikichi

Physics of Plasmas, 24(2), p.022503_1 - 022503_19, 2017/02

https://doi.org/10.1063/1.4975611

The drift kinetic equation describes the collisional (neoclassical) transport in plasmas. Recently, a novel radially-local approximation of the drift kinetic equation, which is called the zero orbit width (ZOW) model, is proposed. In this work, as a numerical verification of the neoclassical transport based on the ZOW model, we perform a series of benchmarks of the neoclassical transport and the parallel flow in three helical magnetic configurations using various types of radially-local approximation models including the ZOW model. We found that the neoclassical transport of the ZOW model can reproduce that based on the other models when the radial electric field and thus the drift is large. Also, it is demonstrated that an unphysical large radial transport, which arises in the neoclassical transport of the other models when the drift is small and compared to the magnetic drift, can be mitigated in the ZOW model.

Radially local approximation of the drift kinetic equation

Sugama, Hideo*; Matsuoka, Seikichi; Satake, Shinsuke*; Kanno, Ryutaro*

Physics of Plasmas, 23(4), p.042502_1 - 042502_11, 2016/04

https://doi.org/10.1063/1.4945618

A novel radially local approximation of the drift kinetic equation is presented. The new drift kinetic equation that includes both and tangential magnetic drift terms is written in the conservative form and it has favorable properties for numerical simulation that any additional terms for particle and energy sources are unnecessary for obtaining stationary solutions under the radially local approximation. These solutions satisfy the intrinsic ambipolarity condition for neoclassical particle fluxes in the presence of quasisymmetry of the magnetic field strength. Also, another radially local drift kinetic equation is presented, from which the positive definiteness of entropy production due to neoclassical transport and Onsager symmetry of neoclassical transport coefficients are derived while it sacrifices the ambipolarity condition for neoclassical particle fluxes in axisymmetric and quasi-symmetric systems.

Analysis of measurements for a Uranium-free core experiment at the BFS-2 critical assembly

Hunter

JNC TN9400 99-049, 74 Pages, 1999/04

JNC-TN9400-99-049.pdf:2.03MB

This document describes a series of calculations that were carried out to model various measurements from the BFS-58-1-I1 experiment. BFS-58-1-I1 was a mock-up of a uranium-free, Pu burning core at BFS-2, a Russian critical assembly operated by IPPE. The experiment measured values of keff, Na void reactivity worth, material sample reactivity worths and reaction rate ratios. The experiments were modelled using a number of different methods. Basic nuclear data was taken from JENDL-3.2, in either 70 or 18 groups. Cross-section data for the various material regions of the assembly were calculated by either SLAROM or CASUP; the heterogeneous structure of the core regions was modelled in these calculations, with 3 different options considered for representing the (essentially 2d) geometry of the assembly components in a 1D cell model. Whole reactor calculations of flux and keff were done using both a diffusion model (CITATION) and a transport model (TWOTRAN2), both using an RZ geometry. Reactivity worths were calculated both directly from differences in keff values and by using the exact perturbation calculations of PERKY and SN-PERT (for CITATION and TWOTRAN2, respectively). Initial calculations included a number of inaccuracies in the assembly representation, a result of communication difficulties between JNC and IPPE; these errors were removed for the final calculations that are presented. Calculations for the experiments have also been carried out in Russia (IPPE) and France (CEA) as part of an international comparison exercise, some of those results are also presented here. The calculated value of keff was 1.1%k/k higher than the measured value, Na void worth C/E values were 1.06; these results were considered to be reasonable. (Discrepancies in certain Na void values were probably due to experimental causes, though the efect should be investigated in any future experiments.) several sample worth values were small compared with calculational uncertaint

Improvement of numerical analysis method for FBR core characteristics (III)

Takeda, Toshikazu*; Kitada, Takanori*; *; *

PNC TJ9605 98-001, 267 Pages, 1998/03

PNC-TJ9605-98-001.pdf:11.65MB

As the improvement of numerical analysis method for FBR core characteristics, studies on several topics have been conducted; multiband method, Monte Carlo perturbation and nodal transport method. This report is composed of the following three parts. Part 1: Improvement of Reaction Rate Calculation Method in the Blanket Region Based on the Multiband Method. A method was developed for precise evaluation of the reaction rate distribution in the blanket region using the multiband method. With the 3-band parameters obtained from the ordinary fitting method, major reaction rates such as U-238 capture, U-235 fission, Pu-239 fission and U-238 fission rate distributions were analyzed. As for the nuclides to be analyzed, the elements of structure material, such as iron, nickel, chrome and sodium were considered. By the present method, all the reactions became larger at the deep region in the blanket. The maximum correction amounted as much as 5%. This tendency lessen the disagreement between the ordinary calculation and the experiment. It was made clear that the treatment in inter-band scattering term is veryimportant because it has large sensitivity on the result. An alternative method to determine the multiband parameters whieh method is based on more direct approach and is free from drawbacks in the present method, was also investigated. Part 2 : Improvement of Estimation Method for Reactivity Based on Monte-Carlo Perturbation Theory. Perturbation theory based on Monte-Carlo perturbation theory have been investigated and introduced into the calculational code. The continuous energy Monte-Carlo perturbation code has been developed by using not only the correlated sampling method which is already used before, but also the derivative operator sampling method. The Monte-Carlo perturbation code was applied to MONJU core and the calculational results were compared to the reference. The change of eigenvalue caused by the change of sodium density in the GEM or dummy ...

Improvement of numerical analysis method for FBR core characteristics (II)

Takeda, Toshikazu*; *; Kitada, Takanori*; *

PNC TJ9605 97-001, 100 Pages, 1997/03

PNC-TJ9605-97-001.pdf:2.82MB

This report is composed of the following two parts and appendix. (I)Improvement of the Method for Evaluating Reactivity Based on Monte Carlo Perturbation Theory (II)Improvement of Nodal Transport Method for 3-D Hexagonal Geometry (Appendix) Effective Cross Section of U Samples for Analyzing Doppler Reactivity in Fast Reactors Part I. Improvement of the Method for Evaluating Reactivity Bascd on Monte Carlo Perturbation Theory. Theoretical formulation in Monte Carlo perturbation method had been checked, and then introduced into a calculation code. The increase of CPU time is about 10 to 20 % compared to that if normal Monte Carlo code, in the cases of same number of history. This Monte Carlo perturbation method found to be effective, because results are almost reasonable and deviations of the results are especially small, by using the Monte Carlo perturbation code. However, there are somc cases that the results of the change of eigenvalues becomes positive or negative by changing the estimator, and there is no reasonable difference in the results between the conventional method, which does not consider the change of neutron source distribution caused by a perturbation, and the new method, which consider that change. Thus it is still necessary to check the Monte Carlo pcrturbation code. Part II. Improvement of Nodal Transport Method for 3-D Hexagonal Geometry It is certain that we can accurately evaluate hexagonal geometry FBR core by nodal transport calculation code for hexagonal-Z geometry named 'NSHEX'. However it is also found that in very heterogeneous core the results is not good enough. Because the treatment of the transverse leakage to the radial direction, which is use for evaluating intra-nodal flux distribution, is not so accurate. For the treatment of the leakage distribution, it is necessary to estimate the nodal vertex flux. In conventional method, the vertex flux estimated by the surrounding node surfacc flux around that vertex. On the contrary,

None

*; *

PNC TJ2222 93-001, 88 Pages, 1993/03

PNC-TJ2222-93-001.pdf:3.54MB

None

Effects of transport theory application on nuclear characteristics of a 1000Mwe LMFBR based on heterogeneous core concept

PNC TJ250 82-13, 44 Pages, 1982/03

PNC-TJ250-82-13.pdf:1.43MB

None

Numerical studies on the neoclassical toroidal viscosity in low collisional plasmas by global kinetic simulations

Matsuoka, Seikichi; Idomura, Yasuhiro; Satake, Shinsuke*

no journal, , 

Effects of non-axisymmetric magnetic field perturbations in tokamak devices in which the magnetic field is axisymmetric have attracted much attention. The non-axisymmetric perturbations produces the so-called neoclassical toroidal viscosity (NTV), which is caused by collisional processes in a plasma. The NTV can make influence on the plasma rotation, which plays a key role in controlling a confinement performance and instabilities. Recently, however, it was pointed out that severe discrepancy exists with regard to the evaluation of the NTV between a Superbanana-plateau theory based on the bounce-averaged particle orbit and a global kinetic simulation. Clarifying physical mechanisms that causes the discrepancy is crucial for precise evaluation/prediction of the NTV. In this work, we perform two different types of global kinetic simulations, of which physical and numerical models are quite different from each other to resolve the issue. As a result, it is demonstrated that the discrepancy arises owing to the lack of the following two mechanisms in the theory; (1) insufficient resonance condition with regard to the precession drift, and (2) the transition processes of the particle orbit caused by the non-axisymmetric perturbations.

Global kinetic simulation study on neoclassical toroidal viscosity

Matsuoka, Seikichi; Idomura, Yasuhiro; Satake, Shinsuke*

no journal, , 

Effects of non-axisymmetric magnetic field perturbations caused by error fields and external perturbations have attracted much attention in many axisymmetric tokamak devices such as JT-60 and ITER from the view point of the control of the plasma confinement performance and instabilities. The non-axisymmetric perturbations cause the neoclassical toroidal viscosity (NTV), and establishing a reliable basis for the evaluation/prediction of the NTV becomes important. Recent studies pointed out that the qualitative discrepancy with regard to the evaluation of the NTV exist between a widely-used superbanana-plateau theory based on the bounce-averaged model and a global kinetic simulation. Hence it is crucial to clarify the cause of the discrepancy. In this work, we perform two different types of global kinetic simulations for the NTV to investigate the discrepancy from the theoretical model. As a result, it is first demonstrated that the discrepancy arises owing to the following two mechanisms in the theory; one is the insufficient resonance condition with regard to the precession drift, and the other is the transition processes of the particle orbit caused by the non-axisymmetric perturbations.